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FOREWORD 


This final report was prepared by the Aerodynamics Section of the Grumman 
Aerospace Corporation, Bethpage, New York for the Flight Mechanics Division, 
Air Force Flight Dynamics Laboratory, Wright-Patterson Air Force Base, Ohio. 
The work was performed under Contract No. F33615-75-C-3073, which was initiated 
under Project No. 1476, "Advanced Wing-Body Aer ¢ Analysis and Design." 
Mr. J. Kenneth Johnson (FXM) was the Project Monitor of this contract. 


The report consists of three volumes. Volume I, entitled "Description of 
Analysis Methods and Applications," describes the methods used to predict sur- 
face pressure distributions and aerodynamic forces on three-dimensional wing- 
body combinations at transonic speeds, including viscous effects. Volume I 
also contains an extensive set of comparisons between numerical predictions 
and experimental results. Detailed instructions required to use the program 
are provided in Volume II, "User's Manual and Code Description." Volume III 
was written at Sybucon, Inc., Atlanta, Georgia and contains a complete descrip- 
tion of the theory and program that computes the full three-dimensional boundary 
layer over the wing. This work was performed by Sybucon under subcontract to 
Grumman Aerospace, Although this program operates independently of the program 
described in Volume II, the input data set required for the full three-dimensional 
boundary layer computation is generated by the code documented in Volume II. 


Mr. F. Berger was the Program Manager; Dr. W. Mason and Mr. D. MacKenzie 
served as Project Engineers. The work was performed in close cooperation with 
the co-authors from the NASA Ames Research Center, Dr. W. F. Ballhaus and 
Ms. J. Frick. Additional contributors to the project included G. Simpers, 

A. Vachris, D. Raila, P. Aidala, M. Sturm and A. Bunnell of Grumman, and 

Drs. F. R. Bailey and T. Holst of NASA Ames. Moreover, contributions have been 
made by A. Chen of Boeing, Drs. R. Melnik, B. Grossman and G. Volpe of the 
Grumman Research Department and Grumman Consultants Prof. A. Jameson, Prof. J. 
Werner and Dr. E. Murman. As noted above, the three-dimensional boundary layer 
program was written by Dr. J. Nash and Dr. R. Scruggs of Sybucon, Inc. 
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This Appendix Volume contains a listing of the program described in 
Volumes I and II. A flowchart of the IBM Version is also included. The 
user should consult Volume I for a description of the method and Volume II 
for the details of program operation and program documentation. 

The flowchart was made automatically by a software package known as 
Autoflow II and available at the Grumman Data Systems. 
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SECTION II 
PROGRAM FLOWCHART 


The Transonic Wing-Body Code has been flowcharted on a Gould Plotter 
using the AUTOFLOW II Program. Because of the size of the program, it could 
not be flowcharted as a single unit. It was broken down into blocks of 
approximately 1000 cards (usually 5 or 6 subroutines). 

AUTOFLOW II generates an internal set of page numbers for each group of 
subroutines. These page numbers are printed in the upper right hand corner 
of each page. Since AUTOFLOW II uses these page numbers to resolve branches 
and subroutine calls (within the same group), they have been left in; if 
they were deleted, a flowchart extending over several pages would become 
virtually useless. 

Branches are denoted by a circle enclosing two numbers; for example, a 
circle enclosing the numbers 15 and 10 would indicate a branch to box 10 on 
page 15 (within the same group). If a decimal number is enclosed in a circle; 
e.g., 16.02, it indicates that the code is continued at box 2 on page 16. A 
branch from another page is indicated by a decimal number; e.g., 17.22, pre- 
ceding an arrow; in this case, the branch would be from box 22 on page 17. 

AUTOFLOW II also produces a listing on non-executable statements (common 
blocks, format statements, etc.) with each subroutine flowcharted. In order 


to conserve space, these pages have been eliminated. 
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7 
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+++++ PRESSURE ORTR 


/FIXCY 168 


BECIN 00 LooP 
UK = KLPIN, 


KL Pit 








WRITE 103 DEY 
¥IA FORMAT 
FROM THE LIST 









eee 





/BLBUTX23/ gi 


BEGIN 0S LooP 
OJ = PSN 












































/BLOUT162/ 

se TO DEV 
103 

FROM THE LIST 


/BLOUT162/ 

LIST = XLENM, 

TLE (Kk), ETA tK) 
CROW 1K) H 






/BLOUTISS/ 
WRITE TO DEV 






/BLOUT 1637 












VIR FORMAT / J Toot vey 
FROM THE LIST / VIR FORNAT / 
FROM THE LIST / 





WRITE 10 DEV 
IBOUT /#BLOUTIBU/ 

¥IA ahs 
t 


ZOCLE = 
0, Sm (ZOCUIK, 1) + 
ZOCLIK, 1)) 


34 





JBLOUTIS?Y 
/ WRITE 14 OEV 


/ 
WIA FORMAT 
101 
FROM THE LIST 
/BL0 i 


LIST = NPTS, NZ, 
KS1A 1 


/BLOUTIS8/ 
WRITE TO DEV 
Fi vie FORMAT f 
/ Sipah neg 
FROM THE LIST /’ 





f WRITE TO DEV 

VIA fGRRAT iy 
104 

FROM THE LIST 

















LIST = 20CLE, 
(Z0CL tk, Dd = 
UL,NPTSHN) 
















/BLOUTISS/ ee 
LIST = NBCO es 
(OITER, LOPAF ves 


ere" AIRFOIL 
ORUINATES /BLOUTX28/ 
va WRITE 10 DEV 


(BOUT ff 
f win FORMAT / 


38 










/EIXCYIG9/ 







BEGIN 00 LOGP f. / 
5 K = KLPIN, / FROM THE LIST / 
KLPOUT 4 





THETA = 
TWIST (K) sRAD 


651 





10/24/77 AUTOFLOW CHRAT SET - FLOCHART BLOUT 


CHART TITLE - SUBROUTINE 61 CUT 


12.39 ro 
/BLOUTIBB/ 
/ WRITE 18 OEY 
f 180UT 
/ ¥1R FORMAT 
f 103 : 
/ FROM THE LIST Of 


BLOUTI6B/ 


LIST = MRCHNO, 
CR PRI TEMP 


AFIUNCTI70/ 
/ WRITE 10 DEY 
/ . PeoUT 
/ S30 FORMAT / 
/ 101 / 
J FROM THE LIST / 


ELC) £20/ 
a 


/BLDUT!70/ 
/ WRITE 10 DEY 
BOUT 


/ WIR FORMAT 
f {01 
Jf FROM THE LUST 


ABLOUTITES | 0 


/ WRITE 10 DEY / 
/ OUT 7 


/ YI FORMAT 
f 103 
¢_fROM THe list / 


ENO OF oo > 


it A Pg 


SEINCTITLS 
BECIN DO LOOP 
5K = KLPIN, /BLOUTI 16/ 
wate 10 OEY 
vIA FeRMaT 
toe 
ee WaITE 10 DEY 
F out 
VIA FORMAT 


/ Loy /FIXCTUSS/ 
Zfron THE LIST 


REWIND TBOUT 


wu FORMATS 


/BLOUTISI/ | 15 
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CHART TITLE - SUBSOUTINE MURMAN IIT, XSHAP) 












/MURMANSS/ 
weamasoe CALCULATE 130 10 
ere eel TON BEHINO 


SHOCK motu neven pew mene me nee s 
ee ee NO 
ENO OF 09 


Lear? 









TS « 
AMAXL Wet K,J,U) 








XS = XSHKPIK,L) - 
*MURFX 287 SHIFT 


(MURFX 364° 
> 


FRLSE” 
tT GT. 







BEGIN 00 Loop 
(60 _K = K7X, KSX 


JO + TENS! - is 


WK. = 0. 





NJ = JTE(K) - 
JLE (KD #2 


NJN = NJ ~ 1 


XS «= 
FRUSHKP IN, LI * 
F)mXOLDIK.L) - 

SHIFT 


#FINCYESS 


XX = (XT (10 - 
i) = x1 40x 


AMURMANI6/ 
FMU(K, JL) EQ. 
FRULK, 






/MURNANIEA S15 


Oo? 


WECIK, JL) = 
NBCLOR K, J) 


- 1b) ¢ 
SySpa adie ax 


/FINCYO61/ 


HOCK AL) 
BCR (K, J) 





"17.01 
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CHART TITLE - SUBROUTINE MURMAN IIT, XSHKP} 


JO= LEM) -1¢ 
Jd 


XX = (XT(J0) - 
XS) /0X 


180 03 


NO 
END OF 00 
LOOP? 


ves (8) 


/MURFX Y3/ [| NOTE 10 


/YANVIS12/ | eae os | 
= 


BED 3 
COEFHTSP BUMP (Xx) 12.14u—— 
NOT 
BEGIN DO LOOP 
30 Ks K2x REX 
190 12 
XOLOIK.LI = 
XSHKP (K, L) 


JAURFX WS¥ 13 
v 


ni NO 


Md 
ENO OF 00 
=§ LOOP? 





JAURFX 417 
COEF = COEU 


PAGE 


17 
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CHART TITLE - SUBROUTINE MUAMAN ITT, XSHKP) 









Fs warre to OEV 
vIA \t wild 
FROM THE LIST 








FX Bs 


laxtes8 OEY 
VIA FORMAT 
LUST 
















/18n_ _73/ | NOTE OS 


LIST + WORD (2ai - 
i ORG (eet) | 





/MURFX Si/ 






WRITE TQ OEV 
VIA FORMRT 






al] aa Ya o 
WRITE 10 DEY 
vIn FORMAT 
FROM re LUST 
















AIBA MOTE 08 





a Va 





Bins 


/MURFX SO/ og 






é 


‘OUT 
UKM, 


p 
“althe 
















af rien 78/ 10 

j WRITE 10 OEY 

vIn Fi 

| VES 

i ‘18K 07 

} ] 

| 1G (26 ARYAN ?/ uv 

; WAITE 10 OEY 
vIAF 

| “cn 






/MURRANTB/ | 18 
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CHART TITLE ~ SUBROUTINE POUT (VAR,NKI,NKO.NJL,NJ2, 1TYPE) } 


coo 
Pour é 
18. 06« 
‘ ou 
NB = (NK2 ~ 
NKLI/10 ¢ 1 
fPour 77 | NOTE 05 


BEGIN 00 LOOP 
WO L +1, NB 









fPOUT 


FALSE.“ 
SSC ne 6. ne > 






f WAITE TO DEV fe 
VIR FORMAT =f 








"e WRITE TO DEV 
Ps VE FeRUnT vA 
“PROM THE LIST / 













/POQUT__18/ | NOTE 16 


LIST © J, 
CYRR UK, IK = ; 
K i 









we on 
SETA MeV SE ad 
WRITE TO Oey / 
vIn Format |” 
; (Ot / 
{FROM THE LIST / 


/ f 

f WRITE TO vey / 

VIR FORMAT / 
108 


/f 








le 










POUT 198" 
e. 


7 
eM OF pO} 
Loar? 






NO 





| 
| 
| 
| sPOUT _1S/ 








/POUT 
L_ne. Le OF as 
~“ ~ 
ee ¢ NOT : 
(1.18) q 
~~ 4 
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Subroutine 


XLDX 
QUIKX 
SAVSOL 
GRAPH 
DSTPLT 


SLETE 


FLOW CHART - 25 


Autoflow Page 
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CHART TITLE - FUNCTION XLDX (KLAST) 


INTEGRATION OF PJUNP 






XLDOX = 
PUMP (KUX, 1) 





/XLDN QU? } 11 


BEGIN Of LooP 
400 K = KROOT 
K™1 






J= JSTECK 


POUR (K) = 
PUMP (K, J) 







18/ 

POUMIKTIP) = 0.0 

EDUM(KTIP) = 
SSPAN 







#XLOX 2D/ 










SINP 

(XLOXX 
EQUA (KROOT) , 
POUR IKROOT! « 
KLAST, JERR} 













#NLOX 21/7 


MLOX = NLDXX 


#XLOX 227 | 08 
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CHART TITLE - SUBROUTINE 


—$$—$_——» 


é quinn / 
—-—— 


6 a SS FROM 
SEMPLE Boot iheur 


SET VALUES OF INDICES 





NREU = JNAXY | 
TBUO = NBS 

NAB = JHAXY 
IBLD = RBS 

NRBS = JAAXK 








AQUIKN 23/ o2 
TBSIQE » LiBuY - } 
RUBLX + 1 | 









EGIN OO LOOP 
OO J = 1, JNAXX 





SRBUID = XINKGD 
YRELL @ XINT ED 
100 


XABS(D = XINX (J) 










5 


Foo 
e e 


Fé 


ETABU IK) = 
= EVAX IN} 


QUIK 
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TRANSFER BO0Y LINE TO 
GRID. aNd GET SLOPE 


/GUTRX 394 


/QUIKX 42 


OATTAN 
(NHBZ, XZWHB, 
we + JHAXX, 


ate 








MPUTE RARER 
Dera TCO 


AF IXELPASY 









A MBITE Yooey / 
/ 
/ via Farmar f 
300 


@FIXELPSS/ 









TQ = 
ETAX (KHODTX) 
208U = 77X (LLBUX) 
ZQ8L = 2TXALBLYO 













JETXELPSS/ 






BEGIN 08 Logr 
= MAXX 


03.28 






B= 2UCLI( - 
2MHBE (17 


A = YMHBI() 


BR = ZMHBL WI - 
ZLCLI tu) 






SAREL = 
QO. 25ePl ele (Bh + 
BS) 











FQUIKX 637 


SRLPHX (Ut = 
SAREL/SBRICK 









ELIPSE 
(Q,ETRBU (KES), 
TQBU, TT. 21 
FNX, PNY, PND) 








JPINELPRU/ 

/ RITE TO DEV / 

VIR_FORMAT 
310 

Fron THE LIST 
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CHART TITLE - SUBROUTINE QUIRKY 


/PINELP on 


, OCBUIKBS! « 
, DISTiV], Ye, 21,22) 
TETRBN (RBS) "= 


ASABE = tty) 
/FUNELPSS / | Qe 

ft ELLIPSE 
(I, ETRBL 2) , 


Z08L,11,21, 
FRX PRY, } 













OCBL UL} 
OISTIL 12,2120) 
/ETRBL (el 






01S a Ae 182) 
tee that RAS! 


a A - 1) 
KBSM1 = ABS - 1 










hoes] Morea 
[isffeit 2 eh] 
| og 


Ae Seen . 


ernst my 172.0 


« fETRBUIKT + 
Fabu ik” = 19) 72,0 


| os aap ws 













alas oe 
ENO OF a 
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OCAUIR? 
OIstur.re ra #2) 
sin = 





= 


iF IXELLOTY \ 


avg AF eS 
7 th : 
NM ane? see 








JFIXELAOB/ | NOTE 1S. 


Ao» (ETABL Ik + 


OCaL IK) = 
Gree etayere te) 


Loor? 


BEGIN O0 LOOP 
ISUK os 2) NBSMI 


+ 
ETABL (KI) 72.0 
= CETABL (kK), ay 
ena tk ~ 1) #20 








ry M15/ iz ; I 


cise 
| 


| Ye, 
72, FNY ie 
_FWE =) 
JFYXEL UI B/ 23 


OCBS (Ll & 


GUSTiVt v2.21, 22) 
7CEEXULLRLY © 
Rie FTN BU 


“PEINELLN 7 | pas 48k 


| Aust 


ax ctBuet, Tiks 
21, PNK PNT 
a FMEI gi! 





(FIXELUISY | 





ah . 
TX (LUBUX < 4)1 
TASMe = LBSTOE - 


PELE THON! 1 NOTE. oe 
Le BEGIN OO Lace | 
Cat 














Bi teneut.t 





Ae (ATX * 1 
ULALN + ETK 
LLALNN fe.0 
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CHART TITLE - SUBROUTINE GUIKY 


ON. 30 
FEIRELI25¢ 


a Seer 


/QUIKX 93/ 


/FIMELIZ7/ | F204 L.) = - 
StL) MENKS 
RECU (20 OCBS tL) ¥FNXS 


K,Jt = - 
DCU (K) »FNXU 


RBCL(20 + KJ) = FFIXE! fs 
DEBL tk) AF AKL : = S 
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CHART TITLE - SUBROUTINE SAVSOL IT TER) 











/SRVSOLIS/ 
JWRITE 10 0€¥ 
! it f 
¥ VIAFORMAT = / 
: 510 
“FROM THE LIST 







J3 oa *Jj- 


XOUT I = 
(XK, J3) ~ 
XLEW (K)) /CORON tK) 


















q 
/SAVSOL2S/ Oo! 

















































1 
/ / sSAVSOL|S/ | NOTE 11 i 
/ REWINO 11 Fi (Uist = «write NOY ig 
f XOUTO(K, J) = END OF DO) > 
i xOUT CS LOOP? 
E /SAVSOL3B/ 12 er a 
/SAVSOL26/ /WiITE 10 OEY TES 
ANRITE = 0 Jrmeggeen f 
FAM THE L(ST 
aes = /SAYSOLS2/ 32 
JSAYSOL2 7/7 03 / car 
/ NAITE TO DEV f 
" VUR_FORMAT —/ 
500 
FROM THE LIST / 









/SAVSOLS2/ 


| LIST = ITER, 
JMAX, KNAX, LMAX, 
TP 


WRITE 10 DEY 












f 


‘i / 
¥IR_FORMAT 
/ "$00 et 
/ FROM THE LIST 








/SAVSOL27/ 


F7GLO = K7X 
__NBQLD = HBX 





J WALTE TO fv / 

/ SIR-FORWAT =f 

ue H ra 
Z_FROM THE LIST / 


/SAVSOLS3/ 34 
1} 













/SRYSOL3B/ 
/ WAITE 10 vey, 
/ 
Pa VIA. FORMAT 
S12 












/SAWSOLUT/ | 
/WAITE TO OEY P 


tt 
ive FORMAT 
Sag 









E} 
LIST © {Xt isis 


























www WALTE / FROM THE LIST of $0) / 1, J ‘ 
Arey p8t 10 SOLUTION Le {FROM THE LIST / fETA(K) RK » 
oar pores 1, KMAN). 
Pree iy Qiu 










L,LWAX) 






ee 


/SAVSOLUT/ 










/SAVSOLII/ 











/ WALTE TO OEV /SAVSOLSUY 36 
VIA FORMAT wo \ 
503 ; TRUE 
FROM THE LIST 
ad 


“nor, goor =) 
Pai 


FALSE {9 
CP) 
8 
JSRVSOLSS/ a 


WAITS ta DEV / 
VIA FeRNat 


soo 7 
é 
0! 
















NX = JTE(K) ~ 
JLE IK) + 1 


NXO KK) = NX 






NX = NXO(K) 












“J 


/SAVSOLSO/ 
WRITE ue DEV 


{ 
VIR FORMAT 

j Sol 
FROM THE LIST 





} 
VIA_ FORMAT 
LIst 















/_ FROM THE ust / 


| 
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CHART TITLE - SUBROUTINE SAVSOL IITER) 











/SAVSOLE7, 19 










UIST ee JNOSE, 
wee hy » 


WRITE 10 DEY 





TRi 





-NOT. BODY 





S01 
FROM THE LIST 












/SAVSOL5S6/ | NOTE G2 
BEGIN 00 LooP 
11 J = JNOSE, 
IT8 


03 


MALTE TO DEV 
VIR WAT 
FROM THE LIST 



















/SAVSOLEB/ 20 







WRITE TO 0 
N = KBMX (J) ‘ in re 
/SAYSOLE: 13 VIA aaa 







FROM THE LIST 







Marre ve ofv 
WIA _FORNAT 
SO! 

FROM THE LIST 












LIST f a ial 


(MBN LO) je 
JNOSEX, JTAILN) 








ENO we S 
Loo! eis 





/SAVSOLS8/ 





*SRYSTLAS/ | NOTE 22 


BEGIN 00 LOOP 
21 J = JNOSEX, 
TRUE AIL 


/SAYSOLE: hy 


















«NOT. EXTASH 









23 19) 5 x2 
FALSE 3 W = KOMXX iJ) Ps 
= © mat. VIS¢ 
is % a 
a 2u 
YES /SAVSOLEY? 15 / WRITE TO DEW |/ 
/ Uh / ie Qh 
WRITE 10 OEY / Vin FORNRT / 
/ 
VIA F t 
pa / farmer / FROM THE LIS 
3 NOTE 07 Z_FROM_THE LIST 
i CONTINUE ] /SaysnL) NOTE 25 errr tT 1 
VISCOUS Sout ion DATA 
/SAYSOLBGs . 6 LIST = Petr tit tT 


{ bur KJ, fem me nem ene 


/SAVSOL60/ o8 LBOXIK, J), 
+ 
WRITE TO OEV = 
1} P} 


vin penn 
LIST 













/SAVSOL? 26 





/SAVSOLES/ 
WRITE 10 DEY 
¥IA F 
FROM 

















etd 


TES 


1 
LIST 
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CHART TITLE - SUBROUTINE SAVSOL IITER) 






WRITE Bi OEY 
¥IA FORMAT 





BEGIN 


/SAVSOL&S/ 
WRITE TO 


2u 
/ WRITE TO EW 

/ 6 vA 

/ ¥1R_FORMRT 


506 Jf 
FROM THE LI st / 


/ WRITE 10 DEV / ssavso.go/_| Ne 










Oo Lor 
QU K = A7K, KBX 


VIA_FORMAT 
S01 
FROM THE LIST _/ 





bev 


16 


TE i? 





Ea] 


en ne 


/ v8 dima / Ir [endif K. 
/ XLS Ne 
ZFOm THe List / Lost 
} 
/Savsmage’ S18 
ENO 4F 00 
~ Loar? 
i ea 
/ Waite 10 ev / 1ES 





J sin sormear/ 
FROM THE LIST 














/SAVSOLES/ 


ons eis s « 


WX = JTE (K) 


JLE(K! + 1 


/SAVSMGI/ | NOTE 13 
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/SAVSOL93/ 
WRITE Ans OEV 
VIA formar 


FROM THE LIST 


/SRYSOL93/ 











WAITE TO OEV f 
VIR_FORMAT 


ta 


va 
a 


Pa 
\. END OF OO CL 
\ Woart 


WAITE TO GEV / 


? Via FORT f 





List = 
(CPU It 


30 









ms DEV 
/ WIR _FORNAT 
Pte v4 
FROM THE LIST 








PRGOE 


10 
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AUTOFLOW CHART SET - FLOCHART SRYSOL PAGE 1! 












Mawcnenm FORMATS 

AMAR 
me Ae me ene ee Se ae ne 
hl kelellalel 





/sawsonot/ 
WRITE 10 DEV 

1 / 
vin, Format 
FROM THE LIST 


Lid? « i 
ACPL th, Bid . 


WRITE 10 DEY 


6 
VIA F 
we / 
Lon ie List / 
























REWIWO If 


FALS! 
OT ati 
7 
"Se 
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CHART TITLE - SUBROUTINE GRAPHICPU,CPL) 


soecenx ASSEMBLE PLOT 
ROUTINE 


ree ee 
PLOTTING AT GRUAMAN 















#GRAOLO20/ 


(8. ,4., 


£GRROL OOS 


sme uheice UPPER AND 
LOWER SU! cr 
DeTnigUTION: AT EACH 


/ GRAD LOSO/ | ou 


CORR = 
AMIN (0. ERROR) 
KLUED = 0 
[STA = 1 
Keel = KTH 


a“ 


/GRAO! 130/ 
KOBL = KROOT 





BEGIN 00 LOOP 
10 K = KROOT, 
KSBL 











Ji = JLE tk} 


WRe Jj - Jie 
XPNT (NXP,K) = 
OOK, J) = 
MLEWIK)) /CORDW (KY 
CPT (NXP) = 
CPUIK, J) 





/GRAO1230/ 14 


CPB INKP) = 
CPL Ik, J} 




















/GRRD1250/ 
CC = CCBAR IK) 
CLL = CLK) /CC 


CHL = 
CMC2 (K) / (CCmCC) 
COL = CO(KI/CC 





/GAAO1290/ \? 


FRAC = 
ETA(K) /YTIP 


(8.5,0.0, - 31 


/SRAO1320 


ENO OF 00 
< CaaP7 


TES 





/GRAD1SYO/ 





PLOT 
(0.,0.,999) 


/GRAD1350/ 


LINENT 
(-1) 


/GRADLIIOS 





PLOT 
10.0,1.5, - 31 


Lc cdetbebal sunrace AND 


Penaath, 


666 
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oNOT. L¥ISC 


SAG Y30, Or ig 
BEGIN DO LOsr ; 
30 K = KROOT, : 

KSBL a 
Ji = JLE(K) 


J2 = JTE (Kt 


/GRAROLYEO/ 


BEGIN 00 Lor 
es Je= Ji, J2 


NXP = J - Sir! 


OSTINXP) = 
OXLSTU (K, NXP) 


OSS (NXP) = 
ORLSTL IK, NXP) 





A ae ES Al rt at ak ee a =e 
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CHART TITLE - SUBROUTINE GRAPH(CPU,CPL) 


wamxnxe FOR USE I 

canes PlotTiNG. IN 
TRIISO serereaen nen rere senercaee 

‘SEE Dene EE ee 


END OF 09 
LOOP? 
/GRRO1T1D/ 
NNX = 31 
OX = NNX 
OX = 1.0/7 (0X - 
3.03 
/GRAROLS6O/ XSTAT (1) = 0. 
XSTAT (2) = 
PLOT G. 25*DX 
(0. .0,.,999) 
/GRAO1760/ 


XSTAT(3) = 
/GAROLS70/ 8. SOxDX 


XSTAT (4) = OX 


/GRROLSIO/ 


PLOT 
0.0.1.5, - 3) 


wane SHIFT 
ZWAX = - 10.0 X-COORDINATES OF 
PLANFORM FOR INPUT TO 
2IN = 10.0 TRILTSO snow nance nse 


WNX = NXP 


BEGIN D0 LOOP 


MPT = MLE[ (1) 


Rees 
WITH EQUAL 
POINTS woe 


(CRAG S30/ 
sein = Fa - 
Yeo! Sse = 
ETA ikRoOT) 
TTREL = YTIP - 
YBASE 


mecca LISE 
KINGH/ARANGE = 
B.O/TTREL. mao acn ann 
I PEE Ope HE DED ene DEE Be ak tim ok 
Lid 

/GRAOLS390/ 

SCALE = YTREL/6. 
XINCH = 
XRANGE / SEALE 
XMOV = XINCH + 
4.0 


YNOV = 
6. /1TREL =I BASE 


/GRAD2030/ 


/GRAD1B80/ XTIC = XMOV - 2.0 
XMAX = XTELCD 


é 
BEGIN DO LOOP 
70_1 = 2, NTEI 


XPT = XTEI(D) 


/GRAN19 10, 


/GRAO20804" 
a 


ENO OF DO 
~ Lorry 


< 


/GRAO1S910/ 
XWAX = XPT 
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CHART TITLE ~ SUBROUTINE GRAPH (CPU,CPL) 


appa ch Carl FOR 
rest ikl Sh ioe 


/GARO2130/ 


/GRRO2160/ 


werrere OEV 
VIA FORMAT 
19 





#GRADZ170/ 
WRITE TO DEV 
VIA FORMAT 
Mt 
09 
REWIND 99 
/GRAOZ210/ 10 







READ FROM DEV 
segs AT 
(NTO THE LIST 





/8RRO2220/ le 


f 
REWIND 98 





/GRRO2Z250/ iY 


S1MBLY 
(-.S, = \ ad 
1,0,.1M, TITLE, 


0.,48) 


/GRROZ270/ 1s 
WRITE 10 DEY 
3g 
VIA FORMAT 
Lie 






/GRROZ290/ 7 

















LIST = (AROD,1 » 
u 


/GAROZ310/ 





/ PEKING 99 
/ 


/GRA02320/ “ee 


SYMBLY 
TMOY ~ 
5, RR, 
4 


uae 27 

FALSE. ee 
Pe 
: 7 





(-.$, - 
1:25, .1 
0. 






wawmwe INTERPOLATE TO 
FIND CP AT NEW POINTS 
FOR GENERAL 

PLANFORN 












/GRAOZWNO/ | 


COUTIJ) = 
CPU IK, J3) 














BEGIN DO LOOP 
165 K = KROOT, 
KTML 





KR = K BF ania + 






TPF (KK) = 
ETA(K) - YBASE 
















7GRAOZNSO/ 


Ad 0 


IAQ2N NOT! 
BEGIN DO LOOP 
INS J 2 it, NNX 


XSTAT(D = 
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/GRADZE1D/ 


END OF OO 
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NO 
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ee oe 


/GRAD2Z710/ 
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t 
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sGRAGZO/ 26 
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XP (2) = (XTET - 
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Z 
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FROM THE LIST / 
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am 
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| PLOT 
2.0, - 4.0, - 
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R 
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YYAL .GT. 
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§. 
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mR 
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| ETA(I) /SCALE 


#SPNOO260/ 
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> 


yk 
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S 
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INTO T fie ey 
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SS 
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ES 
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, 
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ot 


ON 
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A 


a 
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06 
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ak 
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=0 ~ DORAN LINES 
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/ 
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/ 
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a 
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£ 
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2CPPO0220/ 
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1O DEY 
vin ee AT 
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/CPPOOZED/ 
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fe ce 
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VIR FORMAT 
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25, W.R.0. 
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16 
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Srysey 
4.5, .107R.0., 
48) 
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A 
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cee AL FARMS 
CONTOUR PLOT PROGie 


emmame THIS ROUTINE 
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